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Recent decades have seen the discovery of numerous complex materials. At the root 
of the complexity underlying many of these materials lies a large number of possible con- 
tending atomic- and larger-scale configurations and the intricate correlations between their 
constituents. For a detailed understanding, there is a need for tools that enable the detection 
of pertinent structures on all spatial and temporal scales. Towards this end, we suggest a 
new method by invoking ideas from network analysis and information theory. Our method 
efficiently identifies basic unit cells and topological defects in lattices with low disorder and 
may analyze general amorphous systems to identify candidate natural structures where a 
clear definition of order is lacking. This general unbiased detection of physical structure 
does not require a guess as to which of the system properties should be deemed as important 
and may constitute a natural point of departure for further analysis. The method applies to 
both static and dynamic systems. 



I. INTRODUCTION 

Currently, there are no universal tools for examining complex physical systems in a general and 
systematic way that fleshes out their pertinent features from the smallest fundamental unit to the 
largest scale encompassing the entire system. The challenge posed by these complex materials is 
acute and stands in stark contrast to that in simple ordered systems. In crystals, atomic unit 
cells replicate to span the entire system. Historically, the regular shapes of some large-scale single 
crystals were suggested to reflect the existence of an underlying repetitive atomic scale unit cell 
structure long before modern microscopy and the advent of scattering and tunneling techniques. 



'Corresponding author zohar@wuphys.wustl.edu 



2 



This simplicity enables an understanding of many solids in great detail. In complex systems, 
rich new structures may appear on additional intermediate scales. Currently, some of the oldest 
and, after several millennia, heavily investigated complex materials are glasses. Much more recent 
challenges include the high temperature cuprate and pnictide superconductors, heavy fermion 
compounds, and many other compounds including, e.g., the manganites, the vanadates, and the 
ruthenates. These systems exhibit a rich array of behavior including superconductivity and metal 
to insulator transitions, rich magnetic characteristic and incommensurate orders, colossal magneto- 
resistance, orbital orders, and novel transport properties. 

A wealth of experimental and numerical data has been accumulated on such systems. The 
discovery of the salient features in complex materials such as these and more generally of complex 
large scale physical systems across all spatial resolutions may afford clues for a more accurate 
understanding. In disparate arenas, important guesswork needs to be invoked as to which of the 
many features of the physical systems are important and may form the foundation for a detailed 
analysis. With ever-increasing experimental and computational data, such challenges will only 
sharpen in the coming years. There is a need for methods that may pinpoint central features on 
all scales. This works suggests a path towards the solution of this problem in complex amorphous 
materials. A companion work [IJ provides many of the details not provided in this brief summary. 
An explanation of our core idea require a few concepts from the physics of glasses and network 
analysis. Towards that end, we review these concepts below. 



II. A PARADIGM FOR COMPLEX SYSTEMS: GLASSES 



We illustrate the basic premise of our approach by in this work by focusing on glasses. We begin 
by, all too briefly, reviewing a central problem- the detection of natural scales and structures in 
glasses. Complex systems such as glasses are not easy to analyze with conventional theoretical tools. 
As all interactions between the basic constituents of a gas are weak, a gas is easy to understand and 
analyze. At the other extreme, although the interactions in regular periodic solids are generally 
strong, such solids may be characterized by their unit cells and related broken symmetries. The 
situation is, however, radically different for liquids and glasses. Liquids that are rapidly cooled 
("supercooled") below their melting temperature cannot crystallize and instead at sufficiently low 
temperatures become "frozen" on experimental times scales in an amorphous state (a "glass"). 
On supercooling, liquids may veer towards local low energy structures [2, 3j, such as icosahedral 
structures observed in metallic glasses jH [5], before being quenched into this amorphous state. 
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Because of the lack of a simple crystalline reference, the general structure of glasses is notoriously 
difficult to quantify in a meaningful way beyond the very local scales. As such it remains a paradigm 
for the analysis of structure in complex materials. 

The most familiar and oldest technological glasses are the common silicate glasses. More mod- 
ern glasses include phosphate glasses (biomedical applications) , semiconductor chalcogenide glasses 
(optical recording media), and metallic glasses (which have high tensile strength and are stronger 
than steel, resistant to corrosion and wear and tear, have magnetic properties, and an extremely 
high coefficient of restitution). [BHS] Existing work on glasses is vast. Glass formers display several 
key features |TD]. A prominent feature of glass formers is that the viscosity and relaxation times 
can increase by many orders of magnitude over a narrow temperature range. This slowing of dy- 
namics is not accompanied by the thermodynamic signatures of conventional phase transitions nor 
a clearly visible pronounced change of spatial structure. The high number of metastable energy 
states in these systems [HI [12] leads to rich energy landscapes [TlHIZ]. A notable facet of the 
glass transition reflecting structure in the space-time domain is that of "dynamical heterogeneity" 
|18H27j : the fact that the dynamics in supercooled liquids is spatially non-uniform. Many theories 
of glasses, e.g., [10 | I28f[3^ have been advanced over the years. The theory of random first order 
transitions (RFOT) investigates mosaics of local configurations \10\ [28] . As shown in [33j , RFOT 
is related to theories of "locally preferred structures" |33H37j - which, as befits their name, also 
rely on the understanding of natural structures in glasses. Other theories seek a similar quantifica- 
tion of structure. Investigations include spin glass approaches [29j topological defects and kinetic 
constraints [33l ISU EHHlQj , and numerous ingenious approaches summarized in excellent reviews, 
e.g., [141 H2] . There is a proof that a growing static length scale must accompany the diverg- 
ing relaxation times of glass |43]. Some evidence has been found for growing correlation lengths 
(static and those describing dynamic inhomogeneities) |44fl47j . Correlation lengths were studied 
via "point-to-set" correlations [IH] and pattern repetition size [12] • Current common methods of 
characterizing structures that center on an atom or a given link include (a) Voronoi polyhedra, 
[Ml EQI H] , (b) Honeycutt- Andersen indices [52j, and (c) bond orientation [53]. A long-standing 
challenge addressed in this work is the direct detection of structures of general character and scale 
in amorphous physical systems. Towards this end, we briefiy introduce specific concepts from 
network analysis. 
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III. NETWORK ANALYSIS 



Network analysis has been transformative in generating keen new insights in many areas. The 
ideas that we introduce here bring to bear network methods that have been so useful in the social 
and biological sciences to complex physical problems that have not yet been examined before 
through this prism. To address the challenge of detecting and characterizing structure in complex 
systems on all scales, we specifically introduce methods from the new growing physics discipline of 
"community detection" [54J. Our key idea is that any complex physical system may be expressed 
as a network of nodes (e.g., atoms, electrons, etc.) and connecting links. With this representation, 
we may then apply multiresolution methods [55j from network theory to the analysis of complex 
materials. 



Community detection describes the problem of finding clusters- "communities" ~ groups of nodes 
with strong internal connections and weak connections between different clusters (see Figs.([Tl^). 




The definitions of nodes and edges depend on the system being modeled. Between each pair of nodes 
i and j we have an edge weight Vij which may emulate an interaction energy or measured correlation 
between sites i and j [56H58] . The nodes belong to q communities: {Ca}1^i- In our particular 
physics realization, the nodes represent particles and edges their pair-wise interactions. In an ideal 
decomposition of a large graph into completely disjoint communities (groups of particles), there are 
no interactions between different communities; the system is effectively that of an "ideal gas" of 
the decoupled communities of particles. In practice, the task is to find a partition into communities 
which maximally decouple. Such a separation may afford insight into large physical systems. Many 
approaches to community detection exist, e.g., j59H66j and may be invoked in sociology, homeland 
security, and other networks [5 ^1671169] . Two of us earlier developed a method [59] that does not 
suffer from the "resolution limit" that hinders many approaches [70j and that, alongside [61j, has 
been applied to investigate systems with more than a billion links. Earlier works describing the 
method and publicly released code, [SS llTTHTS] contained elaborate and precise technical definitions 
to notions briefly reviewed below. 



A. Partitions of large systems into weakly coupled elements 
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FIG. 1: From [T. A weighted network with 4 natural (strongly connected) communities. The goal in 
community detection is to identify such strongly related clusters of nodes. Solid lines depict weighted 
links corresponding to complimentary or attractive relationships between nodes i and j (denoted by Aij) 
[{Vij—v) < in Eq.([T])]. Gray dashed lines depict missing or repulsive edges (denoted by Bij) [(Vij — v) > 0]. 
In both cases, the relative link weight is indicated by the respective line thicknesses. 

B. Our community detection method in a nutshell 

A parameter 7 in the Hamiltonian below defines the "resolution" of the system. \55\ \7T\ I72j 
We generalize our earlier works by adding a background v and allowing for continuous weights Vij 
instead of discrete weights that are prevalent in graph theory. Our (Potts type) Hamiltonian reads 

^ = 2 E - ^)[^(^ - ^^^^ + - ^)]- (1) 

a=l iJeCa 

In Eq. ([1]), in the inner sum both nodes i and j belong to the same community Ca- The outer 
sum is performed over the q different communities, The number of communities q may be specified 
from the outset or left arbitrary (as in our multi-resolution method below) and have the algorithm 
decide by steadily increasing the number of communities q for which we have low energy solutions. 
\55\ [71] Minimizing this Hamiltonian corresponds to identifying strongly connected clusters of 
nodes. The parameter 7 > tunes the relative weights of the connected and unconnected edges 
and, as advertised above, allows us to vary the targeted scale of the communities sought (the system 
"resolution" ) . As seen from Eq. ([T]) , a high value of 7 leads to forbidding energy penalties unless all 
intra-community nodes "attract" one another [i.e., {Vij — v) < for all (ij) £ Cq], whereas 7 = 
would not penalize the inclusion of any additional nodes in a given community and the lowest energy 
solution generally corresponds to the entire physical system The model for the current application 
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FIG. 2: From [55'. A benchmark demonstrating the accuracy of our method. Left: The benchmark generates 
networks with N = 1000 or 5000 nodes which are assigned to communities of varying sizes specified by a 
power law distribution parameter /?. The fraction of edges that each node has connected to nodes outside its 
own community is fi. The coordination numbers of the nodes are given by a power law distribution specified 
by a. Right: The solutions by our method ("Absolute Potts Model" (APM)) and those found by simulated 
annealing (SA) optimizing a popular cost function ("modularity") are compared to a known answer by the 
Normalized Mutual Information (NMI). Perfect agreement corresponds to a value of NMI = 1. 

could be further generalized by incorporating re-body interactions or correlation functions (such 
as three or four point correlation functions). Details concerning a greedy minimization of Eq. ([T]) 
appear in \55\ \7T\. Somewhat better optimization could be obtained with a heat bath algorithm 
[72] at a cost of a substantially increased computational effort. 

C. Multiresolution network analysis 

We addressed multi-scale partitioning [55] by employing information-theory measures [74H76j 
to examine contending partitions for each system scale. Decreasing 7, we minimize Eq.([T]) with 
progressively lower intra-community edge densities, effectively "zooming out" toward larger struc- 
tures. A key construct in our approach is that of replicas- independent solutions of the same 
problem. This number of replicas p may be set by the user; a higher value of p leads to more 
accurate analysis. In earlier work |55j concerning a static network, replicas were related to one 
another by permuting the numbers of the nodes that form the very same network. However, many 
other definitions of replicas can be considered (as for the analysis of dynamics in complex physical 
systems). We can automatically determine all the natural scales of the system by identifying the 
values of 7 for which these replicas agree most strongly as seen via measures of information theory 
overlaps. We briefly elaborate on these measures in the current context [55J. The probability for 
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a randomly selected node to be in a community a is P{a) = ria/N with the number of nodes 
in community a and the total number of nodes . If there are q communities in a partition 
A, then the Shannon entropy is = —Yla=i 1°§2 ■ The mutual information I{A,B) be- 
tween solutions (partitions) found by two replicas A and B is I {A, B) = Yl^a=i ^6=1 Iv^ ^"''J^ - 
Here, qA and qs are the number of communities in partitions A and B, nab is the number of 
nodes of community a of partition A that are shared with community h of partition B, is 
the number of nodes in community a of partition A, and n^, is the number of nodes in commu- 
nity h of partition B. The variation of information between two partitions A and B is given by 
VI{A, B) = Ha + Hb - 21 (A, B) (0 < VI {A, B) < log2 N). The Normalized Mutual Information 
is defined as NMI{A,B) = l^rfi ^ NMI{A,B) < 1). A high average NMI indicates high 
agreement between different replicas. The VI measures the disparity between different replicas. A 
low VI indicates high agreement between different replicas. A large VI indicates a high variance- 
large fluctuations between the results found in different replicas. The central result of [55j was that 
extrema (including plateaux) of information theory overlaps when averaged over all replica pairs, 
(e.g., the average NMI: 1^ = p(^p__i^ "l^AytB NMI{A, B) and the average VI) indicate the natural 
network scales. |55] That is, we may find the values 7* for which the average Q of information 
theory overlaps such as NMI and VI over all replica pairs [dQ / d'y)\^=Y = and then determine 
the minima of Eq. [Tjfor these 7*(s). The method identifies all "natural" scales of the system. This 
approach is fast [71] and has an accuracy that surpasses methods such as simulated annealing 
(SA) applied to disparate cost functions [551 [7l] (see also Fig. 1). More notably, to our knowledge, 
this approach is the only one that quantitatively evaluates the "natural" partitions over all scales. 
Other current approaches to (non-multi scale) community detection include optimizing "modular- 
ity", |60) dynamics in high dimensions, |64| data compression, [65j. and numerous other ideas. A 
detailed analysis Ref. |59j compared the accuracy of several algorithms for non multi-scale com- 
munity detection. Multiresolution approaches |55 [ I77 | [75] are far more recent. By relatively trivial 
extensions, [1] the method of [55j can be applied to the detection of overlapping communities. 

IV. DETECTION OF MULTI-SCALE STRUCTURES (STATIC AND DYNAMIC) IN 

COMPLEX SYSTEMS 



We wish to analyze complex systems to ascertain general hidden structure in a general manner 
with no prior assumptions as to what the important system properties may be. To achieve this, 
we cast physical systems as networks representing atoms (or electrons etc.) as nodes and setting 



8 





^ 




2 




FIG. 3: From [1 . Left: our replica construction for the physical system at a "static" time to (with no time 
separation between replicas). Right: a similar set of replicas separated by a time At between successive 
replicas. We generate a model network for each replica using the potential energy between the atoms as the 
respective edge weights and then solve each replica independently by minimizing Eq. ([I]) at a given value 
of 7. We then use information measures j55] to evaluate how strongly pairs of replicas agree on the ground 
states of Eq. ([T]) . 

the graph edge weights in Eq.Q to be either (i) pair interaction energies or (ii) experimentally 



measured inter-node (inter-atomic) correlations. As reviewed in Section III C , our approach to 
multi-scale community detection [55j is simple: copies of the community detection problem are 
given to different "solvers" (or "replicas"). If the starting points of different replicas in the com- 
plex energy landscape are different then they will generally arrive at different solutions (different 
community groupings). If many of these solvers strongly agree about some features of the solution, 
then these aspects are more likely to be correct manifesting in extrema of their information theory 
correlations- |55j . When applying this to a physical system, the replicas can be chosen to be copies 
of the system all at the same time in order to detect natural static scales and structures (panel (a) 
of Fig. [3]). Alternatively, the replicas may be copies of the system at different times as in panel 
(b) of Fig. |3] enabling the detection of general spatio-temporal correlations. In both the static and 
dynamic cases, we find the extrema of the information theory correlations as a function of 7 the 
"resolution" parameter in our hamiltonian Eq. ([T]). Once these extremal values of 7 are found, 
the ground states of Eq.([T]) determine the pertinent structures as in [55]. Multiple extrema in 
the information theory correlations suggest multiple relevant length/time scales. In this way, our 
analysis is not limited to the assumption of one or two specific correlation lengths relative to which 
scaling type analysis may be done or what correlation function should be constructed, etc. Rather, 
viable natural scales of the system appear as extrema in the calculation of the direct information 
theory overlaps. 
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V. BENCHMARKS: CRYSTALS, CRYSTALS WITH DEFECTS, AND SPIN SYSTEMS 

Before applying our method to complex systems, we need to make sure that it yields sensible 
results in simple physical cases aside from the standard networks such as those of Fig. [2j Hence, 
we tested it for: (i) Lattices viewed as graphs, (ii) Lattices with defects, and (iii) Defects in spin 
models. 

• (i) Lattices viewed as graphs: Fig. |4] shows our analysis. On the left are the result of our 
multi-resolution analysis. Information theory plateaus correspond to solutions on different scales 
(7). Transitions between different solutions appear as cascades in the information theory measures. 
On the smallest scale (high 7), our approach recognizes the basic units of the lattices. 

• (ii) Lattices with defects: In Fig. [5| we examine a monatomic Lennard-Jones (LJ) system 
(whose ground state is a triangular lattice with inter-particle spacing given by the LJ minimum) 
in which defects in the form of vacancies were inserted. The system is broken into clusters such 
that the defects tend to congregate on boundaries between different clusters. 

• (iii) Defects in spin systems: We investigated Ising systems. In Eq. ([T]), we set Vij to be the 
nearest neighbor bond energies in the Ising model. The multi- resolution analysis led to a cascade 
of structures similar to that in (i) up to the largest domain wall. Partitions into nearly perfect Ising 
domain walls sharply corresponded to VI maxima (and NMI minima). This occurs as the region 
near the domain walls is the one which experiences the largest fluctuations in possible assignment 
to the two bordering domains that it delineates. These benchmark results above are reassuring: 
the systems "find" their natural structures. 

VI. APPLICATIONS TO COMPLEX AMORPHOUS SYSTEMS 

We studied amorphous systems via different approaches starting with that of community de- 
tection In particular, in analyzing complex structures we allowed in some cases for a multiple 
membership of a node in different communities by further replicating each node so as to enable 
overlapping partitions. We used both interaction energies as well as measured pair correlations as 
weights in our analysis. We investigated three different cases: (i) Static and dynamic structures in 
the Kob- Andersen (KA) LJ glass, (ii) Static atomic structures from Reverse Monte Carlo directly 
applied to experimental measurements of Zr8oPt20) and (iii) Dynamic structures in a new model 
amorphous system that emulates experimental results on AlggYyFes. The precise forms of the po- 
tentials in (i) and (iii) are of little importance in the broader context of the method. The goal is 
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FIG. 4: From Multiresolution analysis of a square lattice with periodic boundary conditions treated 
as an unweighted graphs as in earlier analysis |55| : Neighbors have an initial weight of (Vij — u) = 1 in 
Eq. [l]and non- neighbors have an initial weight 0. Left: The panels show the information theoretic overlaps 
between the different replicas when averaged over all replica pairs (see text). These are the variation of 
information (VI), mutual information (I), normalized mutual information (/at), entropy (H) and number of 
clusters (q) in individual partitions. Right: Corresponding partition of the lattice. We use the algorithm 
described at 7 = 60. In this configuration, there were q — 120 clusters with 78 squares, 4 triads, and 38 
dyads which indicates that square configuration dominates the partition, and it shows how our algorithm 
can naturally identify the basic unit cells of the square lattice. 

to analyze amorphous systems and find their natural sub-structures rather than where each atom 
happens to be which constitutes far too detailed and hence useless information. At low tempera- 
tures, we found larger and more pronounced compact structures than those at higher temperature. 
• (i) We studied the glass-forming KA 80:20 binary liquid \80\ [8T] by Molecular Dynamics (MD) 
|82j to simulate an = 2000 atom system. The system is initialized at a temperature T = 5 (in 
the units of \80\ IBlj ) and evolved for a time that is long compared to the caging time. We save 
p high temperature configurations separated by time intervals At of the order of the caging time 
|80j . The system is then rapidly quenched to T = 0.01- well below the glass transition temperature 
of the KA-LJ system. The system is consequently evolved at this lower temperature and again 
we save p configurations separated by the original time interval At. Each of the p copies of the 
system constitutes a replica within the high/low temperature problem. We then employed these 
replicas in our multi-resolution analysis. The results for the low temperature system are shown 
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FIG. 5: From IJ. Left: A plot for the multiresolution analysis of a 2D triangular LJ lattice with periodic 
boundary conditions. The legend for the information theory quantities is as in earlier plots. Edges are 
weighted according the LJ potential, there are two preferred regions, a small peak on the left and a large 
plateau on the right, where the peak here corresponds to the largest possible "natural" clusters. Right: We 
use the algorithm at 7 ~ 31.6 (the left peak) to solve the system. Our method generally places defects near 
the boundaries of the communities in order to minimize the energy cost. 

in Fig. [6} Stable, relatively compact, clusters correspond to the NMI maxima. Additional large 
scale domains in the low temperature system at extrema correspond to a local NMI minimum 
(coincident with a VI maximum) for small 7 ~ 0.2 — 0.3. 

• (ii) Atomic configurations that are consistent with the experimentally determined scattering 
data for liquid Zr8oPt2o [S3HHS] were generated using conventional Reverse Monte Carlo methods 
(RMC) |87H89j . It is notable that KMC analysis of experimental data, unlike MD simulations, is 
not limited by requisite true long equilibration times in MD simulations. A static non-overlapping 
community detection partitioning with the experimentally measured partial pair correlations sub- 
stituting for V in Eq. ([T]) led to extremely large clusters, at a temperature below the liquidus (the 
maximum temperature at which crystals coexist with the melt). The results of the analysis for 
stable clusters are shown in Fig. [7j 

• (iii) Our potential energy functions for AlggYyFes were computed using the techniques of [90] 
employing ab initio results using the Vienna Ab-initio Simulation Package (VASP) |9 11493) . The 
calculated structure factors were compared to experimental data |94j . The potentials are of the 
form V{r) = [(^)°^ + ^ cos [a^r + 04)], with the parameters {aj} depending on the specific types 
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of atom pairs (i.e., Al-Al, Al-Y, ...) (see [T] and, in particular, table 1 therein for tabulated values 
of {tti} for the potentials; further details will be provided in [95j). In Figs. (8|9), we provide the 
structures of Al88Y7Fe5 that our method finds at the two temperatures of T=300 K and T=1500 
K. We found that for a fixed configuration, the community detection Hamiltonian of Eq.Q with 
random edge weights generally exhibits a finite temperature spin-glass type transition [72] (whereas 
for regular lattices Eq.Q is the standard Potts model which exhibits critical or first order transi- 
tions (e.g., a critical transition for g < 4 on a square lattice and a first order transition for q > 4)). 
The structures found are not unique and reflect a conflgurational entropy (different partitions may 
be found for a given value of the resolution parameter 7 that are similar in their overall scale but 
different in precise detail and identities of the nodes). These results may flesh out a facet of the 
glass transition- as the system is supercooled- the effective couplings become quenched. We re- 
iterate that, as noted above, detecting the optimally decoupled structures in more random systems 
such as those that may describe the deeply supercooled liquid constitutes a spin-glass problem. 
|72j The spin-glass transition seen in the community detection problem for a flxed initial atomic 
configuration (as a function of temperature) enhances the change of partitions already evident in 
minimizing the Hamiltonian for the different initial atomic coordinates for the system at different 
temperatures (see, e.g., the larger structures in the low temperature system of Fig. [s] vis a vis 
the smaller structures with much fainter information theory correlations at the higher temperature 
shown in Fig. Mn. 



VII. DISCUSSION 



The detection of structure in general systems is an important problem. We outlined a "first 
principles" network analysis method to ascertain correlations and structures where traditional 
tools of analysis meet with difficulty. With the aid of this method, we were able to detect hard to 
ascertain structures in complex materials such as low temperature glasses. Aside from the canonical 
LJ systems used for simulations of glass former and the more specific (including experimentally 
driven) efforts for the detection of metallic glass systems analyzed here, there are numerous other 
arenas which may profit from such a capability. 

We mention a few of these in passing and speculate on possible applications of our methods 
therein. Omnipresent memory effects and hysteresis appear in amorphous glass systems. [9E\ These 
properties are of direct technological importance for data storage. "Phase-change materials" in the 
form of chalcogenide glasses appear in HD-DVD Blu-ray devices. These devices require materials 
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FIG. 6: From 11]. A three dimensional KA LJ system (see text). Here we allow for overlapping nodes. The 
information overlaps at left show a maximum NMI (/jv) and plateau for other information theory measures 
at (i). Right: We show the typical clusters found. 

that are good glass formers, exhibit large difference in their reflectivity between the amorphous 
and crystalline states, and retain their amorphous state under ambient conditions. The different 
reflectivities of the amorphous and crystalline states allow reading discs by measuring reflected laser 
light. |97H100| Materials similar to Ge8Sb2Teii have these desired properties. Current technological 
progress hinges on vigorous testing of related materials in the hope of eventual improvement. Much 
research has also, in recent years, been done on poly-crystalline Silicon, thin films of poly-crystalline 
CdTe, and similar materials which can be used for photovoltaic cells. [lOlJ A method such as that 
presented here may similarly examine the multi-scale structure of these materials and evolution 
during processing. The basic rudiments of the method can be generally applied to general systems 
having exact or effective interactions between their constituents (liquids, plasmas, etc.). A general 
classical physical system may be represented in terms of a (dynamic) network whose links, loops, 
etc. encode the measured multi-particle correlations and/or interactions. Although we focused in 
this work on real space representations, analysis similar to that outlined in this work may be done 
for nodes that describe the system in Fourier or other spaces. 

We conclude with more speculative remarks about viable extensions of our method to electronic 
and other systems. By attempting to finding an optimal partition into uncorrelated units as we 
have done for the classical systems, a multi-resolution analysis similar to devised here (and in 
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FIG. 7: Left: Information theory measures applied to the Zr8oPt2o system (see text) of = 10000 atoms. 
The computation here were done with V in Eq. [T] replaced by the correlation functions for all different pair 
types (Zr-Zr, Pt-Pt, Zr-Pt) for the RMC data inferred from scattering measurements at 250K below the 
liquidus (1200K). Note the sharp extrema at (i). Right: The corresponding partition. Clusters are assigned 
different colors. 
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FIG. 8: : The result of our community detection analysis applied to Al88Y7Fe5 at a temperature of T = 
300K. The panels at left (a,b) show the information theoretic overlaps between the different replicas when 
averaged over all replica pairs (see text). On the right (c), we highlight the spatial structures corresponding 
to the NMI maximum/ VI minimum. 




FIG. 9: Left: The result of our community detection analysis as applied to Al88Y7Fe5 at a temperature T = 
1500K. When comparing the information theory replica overlaps and structure with Fig. [8j it is evident that 
at higher temperatures, the system is more random. Right: corresponding structures found at information 
theory extrema. 

our more elaborate companion work [IJ) might also be applied to quantum systems where the 
nodes need not conform to localized particles. In the quantum arena, when feasible, a direct 
product representation in terms of decoupled degrees of freedom (and related matrix product states 
|lU2j when these efficiently describe cluster states) may be quite potent. The decomposition into 
optimally decoupled clusters is indeed what the community detection approach seeks to emulate 
in the classical systems that we have examined above. Applications to electronic systems such as 
electronic glasses |103j are natural. Other specific applications may offer insights into the long 
sought diverging length scales (or lack thereof) in strongly correlated electronic systems. In many 
of these systems, the dynamics suggests the existence of a zero temperature quantum critical 
point |104l I105j yet, generally, there are no clear experimental indications of a diverging length 
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scale. Numerous works, e.g., \W(j\ I107j . suggest that a quantum critical point is also present in 
optimally doped high temperature superconductors. Competing orders and multiple low energy 
states may lead to glassy response. [108J As in classical glasses, divergent time scales and various 
non- uniform structures appear in complex electronic materials. These structures are seen by crisp 
scanning tunneling microscope (STM) images and other probes. |109H113] The character of the 
low temperature phases remains a mystery. Generally, a multi-resolution approach similar to 
the one outlined in this work may also be of use for theoretically analyzing general non-uniform 
systems where there are no obvious natural building blocks to consider in performing real space 
renormalization group (RG) type calculations and constructing coarse grained effective theories. 
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